arXiv:1501.07452v2 [physics.optics] 23 Feb 2015 


Purcell factor of Mie resonators featuring electric and magnetic modes 


Xavier Zambrana-Puyalto* and Nicolas Bonod^ 

Aix-Marseille Universite, CNRS, Centrale Marseille, 

Institut Fresnel UMR 7249, 13013 Marseille, France 

We present a modal approach to compute the Purcell factor in Mie resonators exhibiting both 
electric and magnetic resonances. The analytic expressions of the normal modes are used to calculate 
the effective volumes. We show that important features of the effective volume can be predicted 
thanks to the translation-addition coefficients of a displaced dipole. Using our formalism, it is easy 
to see that, in general, the Purcell factor of Mie resonators is not dominated by a single mode, but 
rather by a large superposition. Finally we consider a silicon resonator homogeneously doped with 
electric dipolar emitters, and we show that the average electric Purcell factor dominates over the 
magnetic one. 



INTRODUCTION 


Dielectric Mie resonators made of moderate refractive index materials, such as semi-conductor materials in 
the optical or near infra-red regime, can feature electric and magnetic modes of low order. These modes can 
be easily and efficiently excited either from near or far field illumination [1 -3]. This novel class of photonic 
resonators is very promising to design directive antennas [4-6], to enhance the electric or magnetic near field 
intensities [7-10], to design subwavelength sized light cavities [8, 11, 12], to host frequency mixing processes 
such as third harmonic generation [13], and even to create isotropically polarized speckle patterns [14]. 
Contrarily to spherical nano-metallic particles, dielectric Mie resonators exhibit both electric and magnetic 
resonant modes [2, 3, 15]. However, it is still unclear whether magnetic modes are a good platform to enhance 
light-matter interaction in silicon subwavelength-sized cavities. An interesting way to address this question 
is to study the Purcell factor of these cavities. 

The Purcell factor is a figure of merit widely used to characterize light-matter interaction in photonic 
cavities [16]. In quantum electrodynamics, for example, it has been widely used to study Rabi oscillations 
in the weak-coupling regime [17]. The Purcell factor quantifies the enhancement of the spontaneous decay 
rate of a dipolar emitter coupled with a cavity, and it is defined as [16, 17]: 



with ui being the frequency of the dipole; c the speed of light; Q M the quality factor of the electromagnetic 
mode of the cavity (E^) that the emitter is mainly coupled to; U M the effective volume of the mode E M ; and 
r (ro) the spontaneous decay rate of a dipolar emitter in the presence of an optical cavity (vacuum). 

The design of photonic cavities at the nanoscale brought the attention of the nano-optics community to 
the Purcell factor [18]. Nevertheless, the first attempts of using the Purcell factor given by eq.(l) to compute 
the enhancement of spontaneous decay rates (P/To) were not as accurate as expected. It was seen that the 
Dyadic Green tensor formulation [19] used to compute P/To did not yield the same results as the Purcell 
factor given by eq.(l). As noted by Koenderink [20], the use of eq.(l) to compute spontaneous decay rates 
is based on a few restrictive assumptions, which are not usually fulfilled by photonic nano-cavities. Hence, 
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some changes needed to be done in the Purcell factor formula so that Fp(u>) = T/To- These changes can 
be summarized in two: extension of the Purcell factor to a scenario where the dipolar emitter couples to a 
superposition of modes of the cavity; and redefinition of the effective volume V fl of each of the normal 
modes of the resonator. One of the first attempts to redefine the effective volume for spherical cavities was 
proposed by Colas des Francs et al. in 2012 [21, 22]. The proposed method used the equality between the 
classical normalized decay rates [23-25] and the Purcell factor given by eq.(l). That is, the effective volume 
Vfj, was computed by substitution: 


6t rc 3 ^ T 6nc 3 Q^ T 0 

Fp{u) = = K => ^ = 


( 2 ) 


With this empirical formulation, an effective volume is obtained for each multipolar order of the spherical 
resonator. Recently, attention was focused on the definition of the effective volume with respect to the normal 
modes of arbitrary resonators [26-28]. The normal modes of the system are solutions of Maxwell equations in 
the absence of sources, with outgoing radiation conditions, and with a certain normalization [26, 28-30]. Due 
to the outgoing radiation conditions, the Hamiltonian of the system is non-hermitian and the normal modes 
of the system have complex eigenfrequencies [26, 31, 32]. Let us notice that a normalization condition 
of the kind f v |E M (r, u„)|W cannot be applied to normal modes with resonant complex frequencies 
because E M (|r| —> oo—> oo [20, 26-29, 33]. In this work, we define the normal modes associated to the 
eigenfrequencies following the normalization condition given in [28, 34]: 


1 = / ’ e ( r ) E A‘ clr + 7p-y 


JdV 


„ 9 9E^ 

E„ • — r 


— r 


9E m 3E m 


AS 


( 3 ) 


^ ds dr ' dr ds 

where all the magnitudes that depend on the frequency are particularized at w = At this point, two 


comments are in order. Firstly, the normalization condition given by eq.(3) can also be expressed as [26, 28]: 


1= f [E m • e(r)E„ - • Mo H m ] dr (4) 

Jv 

As Sauvan et al. showed in [26], this type of normalization is especially well-suited for numerical calculations. 
The underlying reason is that the volume integrals can be numerically evaluated using a perfectly-matched- 
layer technique. Furthermore, their normalization condition is easily extended for dispersive and magnetic 
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materials [26]: 


1 = 


d{ue)^ d{ujn) u ' 


dr 


(5) 


Secondly, the basis of normal modes (E M ,H M ) for leaky cavities with arbitrary geometries is not complete, 
in general [27]. However, for spherical leaky cavities, it has been demonstrated that the multipolar fields are 
a complete basis to describe the problem [35, 36]. 

GENERAL EXPRESSIONS FOR ELECTRIC AND MAGNETIC EMITTERS 

The recent advances on the modal analysis of Purcell factors were performed with electric dipolar sources. 
The formulation has been presented in [26-28], where the enhancement of the decay rate was computed as: 


F(u) = —Imjp* • t? (e) (r 0 ,r 0 ) • p} 
nuj V ) 


37rc 3 




1 


- uj) 


( 6 ) 


with 


y,= 


(7) 


(u p ■ E jU (r p , w^)) 

where p = pu p is the electric dipole moment of the emitter, and u p is a unitary vector. However, interest in 
magnetic spontaneous emission has been growing over the last years, in particular using trivalent lanthanide 
ions [37-39]. This strong interest in coupling lanthanide ions to nano-cavities makes us emphasize the fact 
that both electric and magnetic emitters can give rise to an enhanced spontaneous decay rate. Following 
the formalism shown in [28], next we give the Purcell factor formula for an emitting magnetic dipole. The 
interaction Hamiltonian for a dipolar magnetic interaction is [40]: 


Hr = —m • B 


( 8 ) 


where m = ?nu m is the magnetic dipole moment of the emitter, and u m is a unitary vector. If is the 

magnetic dyadic Green function [41], then the Purcell factor can be computed as: 


F(u>) = -^4m (m* • (r 0 , r 0 ) • ml 

nui y J 


37rc 3 


J2 lm 


VpWniUn - w) 


(9) 
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where 




(u m • B M (r p ,w A ))" 


(10) 


It is important to note that eqs.(6) and (9) cannot be directly added. That is, the enhancement of decay 
rates of a dipolar emitter which is a linear combination of an electric and a magnetic dipole (e = ap + /3m) 
is computed with the following formula [19]: 


/ Re{f-E}dy (11) 

M) zl 0 Jv 

where j = — iojap5(r — r p ) + V x /3m <5(r — r p ) [42, 43], and P 0 is given in [19]. Using the currents produced 
by both dipoles, P can be written as: 


P = |lm{a V • E(r„) + /3*m* • B(r 0 )} (12) 

where E(ro) and B(ro) can be computed using the Green tensor formalism [19]. Notice that when the 
emitter is a superposition of an electric and a magnetic dipole, E(r 0 ) ^ G ( e )(r 0 , r 0 ) • p since: 

E(r 0 ) = iu>nno [ t? (e) (r, r 0 ) ■ j(r)dU (13) 

Jv 

with j(r) given by the equation above. 


MULTIPOLAR EXPRESSION OF THE NORMAL MODES 

Here, thanks to the analytic formulation of Mie theory, we analytically compute the normal modes of the 
system and their associated normal frequencies. In particular, for a Mie resonator, due to the symmetries of 
the system [44], the normal modes of the system are the multipolar fields, E M = E, where j = 1, ..,oo, 
m z = q = l,..,oo and (x) = (e),(m). The meaning of j,m z , and (x) is directly related to the 
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following differential operators [44-46]: 


J 2 

Jz 

n 


(x) 

j,m z ,q 

(x) 

j,m z ,q 


= JU + 1)E 


(x) 

j,m z ,q 


(x) 


= m ^j,rn z ,q 


(-)j E (m) 


(_)t+ 1 E (e) 


(14) 

(15) 

(16) 


J,m z ,q 


Nonetheless, q is different in this sense. Given a fixed combination of (a/)}, there exists a numerable 

(x') 

infinite of q values such that E), ^ is a normal mode of Maxwell equations with an associated eigenfre- 
quency = w), q [21, 47]. These eigenfrequencies w), q can be found by solving the following two 
transcendental equations: 



h f ] (“S*. ,q) K*£L , q 3i ,<,)]' 

for 

(x) = (e) 

(17) 


h f ] ( X j^ z ,q)l n r X jH z , q jj Kz-i,,)]' 

for 

(x) = (to) 

(18) 

i ( x ) R ( x ) 

where x)^ z q = 






(x) (x) 

Here, it is important to note a few facts. Firstly, eqs.(17-18) do not depend on m z , therefore = uj]j q . 

Secondly, in eqs. (17-18), the dependence on j is explicit on the spherical Bessel functions of the different 
kinds, and the dependence on (x) is apparent as the equations for electric and magnetic frequencies are 
different. However, the dependence on q is not explicit. That is because, given a fixed j', q is an integer that 
numerates all the solutions to eqs.(17-18) in ascending order, i.e. Re j < Re |w jf j < Re jtujfg j < .... 

Due to the normalization condition imposed by eq.(3), the normal modes of the system are [28, 34]: 


E 


(m) 


J (r,0,<j>) = A^R j (k^) 


1 dY jtm 

sin 0 dcj) 


dY~ 


dO 


A, 


= “T - 7777“ \j(j + m(^)Y j , mz r+ 


e(r)kj^r L 

'^R^q) 


]\ n ‘j,q) ± 3,‘ 


dYi 


3,m z 




1 dY 


J,m z 


dr 


dO 


dr 


sin 0 d(f> 


(19) 

( 20 ) 


where (r, 0, <p) are the spherical coordinates in the real space; k ) T ™- ) and kffl are related to the eigenfrequencies 
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(x) (x) 

of the mode, k) = ud ’/c with c being the speed of light in vacuum; 


-4<”> = 


j(j + l)f? 3 («y - 1) 


A $ e) = nlA^ 


\ 


jj-iinrkj^R) j 


jjinrk^R) n r k^ q R\ <h {e) m2 


Wi 


+ 


j(j + 1) 


(*&) 2 


are normalization constants, with jj(x) being the spherical Bessel function; 

jj{n r k^r) 


Ri& = { 


jj( n rk^qR) 

hj(k%R) 


for r < R 


for r > R 


where hj(x) is a spherical Hankel function of the first kind; 


Y = 


/2 i + Ki-KI )\ P ^\ {cose)XmM 


2 (j + \ m z\V- 0 

are the real spherical harmonics, with p| m ’(cos 0) being the associated Legendre functions, and 

sin (m z <fi) 


Xm z — \ 


V 77 

1 

cos (m z 0) 


for m < 0 
for m = 0 
for m > 0 


being an azimuthal function; finally, the electric permittivity function is 


-O’) = 


ii for r < R 


(21) 

( 22 ) 


(23) 


(24) 


(25) 


(26) 


1 for r > R 

Notice that a constant magnetic permeability /i = 1 is assumed both for the embedding medium and the 
sphere. 


Analysis of the effective volume 

Here, we show that the effective volumes of Mie resonators can be related to the translation-addition 
coefficients , 4'Vm' } an< ^ t ^ ie coefficients {aj>,bj>} (see Appendix). First of all, we re-write the Purcell 
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factor formula given by eq.(6) using the normal modes defined in eqs. (19-26) and the eigenfrequencies given 
by eqs.(17-18): 


= — £ 


LU 


I I II • 


v ie) J e 


e) ( (e 

,1 


e) 


Im 


( m) f 
3,1 


V> ' LU' ' I LU' ' — LU 
3,m z ,q^j,q \^J,q w 


(27) 


where 


V, 


(*) 


3,m z ,q 


(28) 


(u p . E w,yr p .^-')j 

That is, ten parameters need to be specified to compute the effective volume for a dipolar emitter using 
eq.(28): = ^ <ig (r pi Up). 

Now, it is enlightening to compare eq.(27) with the classical decay rates of the same system, which can be 
found in [25]: 

p 3 c 6 

H [\ a E^,m z ) + aja%(j,m z )\ 2 + \a d M {j,m z ) + b ja d M (j,m z )\ 2 ] (29) 


P( 


J,rn z 


The multipolar coefficients a%(j,m z ) and m z ) defined in [25] modulate the dipolar field exactly in the 
same manner as the translation-addition coefficients do in eq.(A.l) [25]. They are defined as: 


tj% z = oc Up ■ A^(r p ,w) 

4^1 = a d M (j,m z ) oc Up ■ A^(r p ,w) 


(30) 

(31) 


Hence, by comparison, we observe that the effective volumes can be analyzed thanks to the translation- 
addition coefficients: 


\tj%, m+ a 3 m i 2 <— Im ' 


1 


(32) 


T/ (e) , ,(e) / (e) _ \ 

1 ( ''j,m z ,q U, j,q 

It is observed that the Purcell factor dependence on the azimuthal number m z is fully captured by the 
effective volume. For instance, if (w) = 0 => Vj^, — > oo. But the dependence on j is more complex, 
and so is the dependence on to, which is given by a summation of the kind l/(A g — B q ui). Thus, the 
effective volume quantifies the coupling between the electromagnetic field created by the emitter and the 
normal mode in consideration, as well as their spectral overlap. 



PURCELL FACTOR OF SILICON NANO-CAVITIES 

For the rest of the study, we consider a silicon-made Mie resonator with a refractive index n r = 3.5, 
magnetic permability /r = 1, of radius R = 0.2 fj,m, embedded in air. The normal modes of the system and 
their eigenfrequencies are computed using the formulas given by eqs.(17-26). It will be shown that even 
in the more symmetric situations, the Purcell factor has contributions from more than one normal mode. 
Also, we will show that, thanks to symmetry considerations, it is possible to predict the value of for 

some cases. Finally, we will show that given a multipolar order j', the average electric effective volume (it 
is defined next) is always smaller than the average magnetic effective volume. 

Dipole in the center of the Mie resonator 

An electric dipole is first supposed to be placed at the origin of the resonator, oriented along the z axis, 
i.e. p = pi. Let us remind that the emitting frequency does not need to be given to compute the effective 
volumes. Thus, the electromagnetic field emitted by this dipole can be expressed simply as E a = A) q, i.e. 
a single multipole (see appendix). That implies: 


tj£ z = 0 , V (j, m z ) 

= 0 , V (j, m z ) ^ (1,0) 

Then, using eq.(32), it is observed that 

v j^l, q 00 - V(j,m 2 ,g) 

^ 00 > V 0>U,g) ^ (1,0, g) 

Thus, the Purcell factor is given by: 


(33) 


(34) 


37TC 3 

F(u) = —2^ lm 

9=1 


V (e) (e) ( (e) _ \ 

' / l, 0 ,q U} l,q \^ W 1,9 


In Table I, we give the first four effective volumes and frequencies corresponding to q = 1,.., 4 units of c = 1. 


Using these values, the Purcell factor F(u>) as a function of the frequency of the emitter (w) is displayed in 
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TABLE I. Effective volumes and eigenfrequencies associated to modes of the kind Ej®q , for q = 1,4. The effective 
volumes {V}*q q ) are computed with respect to an electric dipole oriented along the z axis and located in the center of 
the Mie resonator, using eq.(28). The resonator is embedded in air, has an index of refraction n r = 3.5, and a radius 
R = 0.2pm. The eigenfrequencies (w^) associated to are computed with eq.(17). 



<7 = 1 

q = 2 

q = 3 

<7 = 4 

y( e ) 

V l,0 ,q 

0.2514 + 0.2494* 

0.0383 + 0.0111* 

0.0151 + 0.0011* 

0.0077 + 0.0004* 

w i,« 

5.2940 - 2.4525* 

6.0470 - 0.6070* 

10.9925 - 0.5150* 

15.5635 - 0.4665* 


Fig.(l). As it can be observed, F(uj) has contributions from different modes. When one of them is resonant, 
the behavior can almost be described with only one mode, but when they are not, this is no longer true. In 
addition, as it was shown in [20], when the Q-factor of a mode is not very high (as it is the case for E^ e g -J, 
the Purcell factor can be underestimated if only one mode is considered. 

Dipole displaced along the z axis 

Here, we consider the case where an electric dipole with p = pz is translated from the center of the 
resonator to its surface along the 2 axis. As discussed before, the electric field emitted by the dipole from a 
frame of reference centered on itself is given by E d = A^q. However, the description of the normal modes 
of the cavity with respect to a frame of reference centered at the displaced dipole becomes cumbersome. 
It is easier to describe E d from the frame of reference centered at the center of the cavity. Therefore, the 
translation-addition coefficients are used to express E d with respect to a different point Zq on the z axis. 
As we showed in eq.(A.l), the expression is a general superposition of multipolar fields. Now, we demonstrate 
that a translation along the z axis restricts this general superposition to a case where = f(j)^o,m z and 
therefore: 

Ed = A i?o( r - 2 oz) = 4 ',o A t,o ( 36 ) 

j' 

On one hand, the multipolar fields are eigenstates of the z component of the angular momentum, J z , i.e. 

= mAjl [44, 46, 49]. On the other hand, a translation along the z axis can be expressed as a 
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Purcell Factor as a function of q and co 



FIG. 1. Purcell factor (F(to)) of an electric dipole oriented along the 2 axis located in the center of the resonator. 
Both the total Purcell factor and its partial ^-contributions shown by eq.(35) are depicted. 


function of the 2 component of the linear momentum, P z : T z (Az) = exp (— iP z Az). Then, because J z and 
P z commute, [J Z ,P Z ] = 0, a translation along the 2 axis must maintain the 2 component of the angular 
momentum [46, 50]. A justification why t^l = 0 can be found in [46]. With the relations brought forward 
by eq.(32), we obtain that: 


V, 


(m) 


V, 


J,m z ,q 

(e) 

j,m z ,q 


OO , V(j,m z ,q) 

00 , y(j,m z ,q) ^ (j, 0, g) 


(37) 


Hence, the Purcell factor can be analytically written in this occasion as: 




T/( e ) 11 ^ 
V j,0,q U} j,q 




(38) 


Notice that even though the formula is still very simplified (with respect to eq.(27)), the number of mode 
volumes and eigenfrequencies that need to be computed in order to calculate F(u) increases a lot. In Table 
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TABLE II. Eigenfrequencies of a Mie resonator embedded in air, with n r = 3.5, and R = 0.2/rm. Only the frequencies 


in interval Re j € [4,16] are listed. 



q = i 

q = 

2 

q = 

3 

q = 

4 

, ,( e ) 

W l,9 

5.2940 - 2.45251 

6.0470 - 

0.6070* 

10.9925 - 

0.5150* 

15.5635 - 

0.4665* 

, .(e) 
w 2 ,9 

7.7377- 0.1466i 

10.5360 - 

3.3610* 

12.7485 - 

0.6925* 



, .(e) 

U 3,q 

9.6137 — 0.029751 

14.2235 - 

0.4400* 

15.8325 - 

3.8870* 



, ,( e ) 

W 4, 9 

11.392145 — 0.00597951 

16.1388499 - 

-0.138451* 





,, (e) 

W 5, 9 

13.1071675 — 0.00119451 







. .(e) 

W 6,9 

14.78344639 - 0.0002351* 








II, we give the eigenfrequencies of all the normal modes whose frequencies lay in the interval given by the 
frequencies given in Table I, i.e. Re j € [4,16]. In Fig.2, we plot log (F(w, £o))> *-e. the Purcell factor 
as a function of the emitting frequency of the dipole (w) and its position in the z axis (zq). The logarithmic 
function is used to visualize the behavior of log (F(lj,zq)). 


It is observed that when the dipole is very close to the origin, log (F(u, Zq)) has broad resonances, as 
obtained in Fig.l. However, when the dipole is displaced from the origin, its coupling to the resonator gets 
stronger as well as more resonant. Note that even though some of the resonances have an extremely narrow 
linewidth, their contribution to the Purcell factor is not present for all the positions of the dipole. For 
example, almost does not contribute to the Purcell factor for ^ < 0.3. The reason why this happens is 
that the effective volume V§ e 0 \ is very large for ^ > 0.3. This can be observed in Fig.3, where the effective 
volumes of all the modes used to compute F(u>,z o) in Fig.2 are given. Finally, note that even though the 
most resonant behavior is given by the modes with q = 1, the whole Purcell factor picture would be different 
if only modes with q = 1 were included. A comparison between the two has been depicted in Fig.4. 
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to (rad/s) 


FIG. 2. Logarithm of the Purcell factor (log [F(lo, zo)]) as a function of the frequency u> of the emitter and its position 
zo in the z axis. The orientation of the dipole is along the z axis. The Mie resonator is embedded in air and is defined 
by n r = 3.5, and R = 0.2 fim. 


Effective volumes. Average values 


In this section, we describe the computation of average effective volumes and show that given a multipolar 
order j, the average effective volumes are always smaller for electric modes. Suppose that a Mie resonator 
is doped with a large number of dipolar emitters. In general, both their locations and orientations in the 
resonator are unknown. Thus, we can define the average effective volume as: 


y.( x) 

3,m z ,q 


3r 2 d?’dfl 


In [r • E • l 2 ■ 


[0-E 




[0- E 


(*) 12 
j,m z ,ql 


(39) 


Because the order of the integrals can be interchanged, we can also define an average angular effective volume: 


V, 




r)) = 


3dfi 


In ff ■ Ey 


x) 


(r)] 2 + [fl-E^ i9 (r)] 2 + [0.E^ Mi(| (r)] 


(x) 


(40) 
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Effetive volumes as a function of z Q 




R 


FIG. 3. Logarithm of the effective volumes (log [Vj'oqj) of normal modes of the kind Ej e Q q as a function of the 
position of the emitter z o- Only the normal modes associated to the eigenfrequencies given by Table II are considered. 
The emitter is an electric dipole oriented along the 2 axis, and its position is moved along the z axis. The Mie resonator 
is embedded in air and is defined by n r = 3.5, and R = 0.2 fim. 


so that the average effective volume is computed as: 


( V tm z , g ) = j o r 2 (vg^(r))dr 


(41) 


(’a;') (x) 

Note that in this case, there is only one symmetry argument, which is that Yj„ lz q = . The reason 

stems from the definition of the normal modes of the systems given by eqs.(19-26). In particular, the only 
dependence of El _ in m z comes via Y) m , and Xm z i an d it can t> e seen that an integral over 27 t equalizes 
the behavior of the modes with m z and —m 2 . Thus, in general, the average effective volume (vj^ q ^ 
does not tend to infinity for all j,m z ,q, ( x ). In Table III we present the values of the inverse of 
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FIG. 4. Comparison between a) the Purcell factor computed using all the normal q-modes of the resonator and b) 
the Purcell factor computed only considering modes with q = 1. In both cases, the Purcell factor is computed as a 
function of the frequency uj of the emitter and its position zo in the 2 axis. The emitter is oriented along the 2 axis 
in both cases. The Mie resonator is consideration is embedded in air and is defined by n r = 3.5, and R = 0.2 fim. 

for different electric and magnetic modes with q = 1. If 1/ q ^ is large, the mode in consideration 

contributes significantly to the Purcell factor at w = Rc j. That is, a large value of 1/ (v^ z q ^ means 
good coupling between that mode and the resonator. Two observations can be made. Firstly, it is seen that 
the coupling 1/ q ^ decreases with the multipolar order j. This result can be interpreted using the 

radial dependence of multipolar fields. Indeed, the intensity of a multipolar field has a radial dependence 
that is given by a spherical Bessel function jj{kr). Now, when given a fixed kr* position close to the origin, 
jj-i{kr*) > jj(kr*) [51]. That is, lower multipolar orders have a non-null intensity in the center of the 
sphere, whereas the intensity of higher orders is pushed to the surface. Consequently, smaller values of 
I- <jl contribute to the integral given by eq.(39) for most of the points in the domain. 

Secondly, given a multipolar order j, the coupling of the magnetic modes is always smaller than the coupling 
to the electric ones. This fact can be understood using eq.(32). As we saw for the case of a displaced dipole 
along the 2 axis, the coupling of a dipolar emitter to magnetic modes is completely suppressed. This fact 
changes when the dipoles are not located or oriented along the z axis, but the coupling with electric modes 
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TABLE III. Computation of 1 / (vjm z q ^ f° r different normal modes of the resonator. All the volumes have been 
computed with q = 1. The variation in {j, m z , (x)} is made explicit in the table. Large values of 11 / 
indicate good coupling between the uniform distribution of dipoles in the resonator and the mode in consideration. 
The resonator is embedded in air and is defined by n r = 3.5, and R = 0.2 \xm. 


j 

m z 

(e) 

(m) 


m z = 0 

22.6754 - 12.5357* 

0.778488 - 0.0930585* 

j -1 

m z = 1 

22.6754 - 12.5359* 

4.69505 + -0.561234* 


m z — 0 

8.86982 - 0.944547* 

0.0222563 + -0.0020265* 

3= 2 

m z = 1 

2.89243 - 0.745228* 

0.310968 - 0.0283146* 


m z — 2 

0.674872 +-1.2691* 

0.0814357 - 0.00741497* 


m z = 0 

4.39869 - 0.0542868* 

2.14545 x 10- 5 - 6.63047 x 10 -7 * 


m z — 1 

0.0743161 + 0.120752* 

0.000216611 - 6.69432 x 10" 6 * 

j = 3 

m z = 2 

-0.157549 + 0.0998519* 

0.000122304 - 3.7798 x 10~ 6 * 


m z — 3 

-8.03716 x 10~ 10 + 2.58389 x 10" 10 * 

1.60496 x 10 -12 - 4.96009 x 10 -14 * 


(given by the translation-addition coefficients t ) still remains stronger. 

Finally, we plot 1/ (Vj^ q {r )^> for a few normal modes. This is depicted in Fig.5 and Fig.6. Note that 
the results of Table III are computed as radial integrals of these plots. 


CONCLUSIONS 

In conclusion, we derive the Purcell factor for magnetic and electric dipolar emitters. We particularize it 
for the case of spherical Mie resonators and emitting electric dipoles. Using the eigenfrequencies and the 
analytic expressions of the normal modes of the system, we analyze the meaning of the effective volume 
for spherical cavities. It is seen that its value is related to the multipolar coefficients in the translation- 
addition formulas. That is, the effective volume of a normal mode locally quantifies the strength of the 
coupling of a dipolar emitter to that normal mode. Finally, the Purcell factor is computed for different 
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(j-m z )=(l,0) (j,m z )=(l,l) 


Effective volume, electric modes Effective volume, electric modes 



FIG. 5. Absolute values of the inverse of the angular average effective volume (| 1/ q (ryj \ ) for j = 1 and q = 1. 

Large values of 11 / (v^_ ? (r)^ | indicate good coupling between the uniform distribution of dipoles in the resonator 
and the mode in consideration. The resonator is embedded in air and is defined by n r = 3.5, and R = 0.2/rm. 

symmetric scenarios. We observe that the magnetic effective volumes are —> oo for some symmetric 

configurations, showing no coupling between the electric emitter and the magnetic modes of the structure. 
However, in less symmetric configurations, this is not necessarily the case. Finally, an average effective 
volume is computed both for electric and magnetic modes, revealing that i.e. the 

coupling to electric modes is generally better than the coupling to magnetic modes. The situation could be 
reverted if the effective volumes were computed with dipolar magnetic emitters. 
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(j,m z )=(2,0) (j,m z )=(2,l) (j,m z )=(2,2) 



FIG. 6. Absolute value of the inverse of the angular average effective volume (| 1/ q (r | ) for j = 2 and q = 1. 

Large values of 11 / (v^ 5 ( r )^) I indicate good coupling between the uniform distribution of dipoles in the resonator 
and the mode in consideration. The resonator is embedded in air and is defined by n r = 3.5, and R = 0.2/im. 
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APPENDIX 


Here, we give describe the interaction of a dipolar emitter and a Mie resonator. As it is described in [52], 
an emitting dipole is a multipolar field with j = 1. If the dipole is electric, then its associated electric field 
can be described as A^ with respect to a frame of reference centred on the dipole itself, where m z depends 
on the orientation of the dipole. In fact, m z = ±1 if the electric dipole is given by p = , or m z = 0 if 

p = pz. Then, when Aj m is expressed with respect to a reference of frame which at a distance ro of it, the 
dipolar field is re-expressed as 


E d = A (e) 


(V-rvA— A e > A( e ) i +'l m ) A ( m ) 

l,m z v I 1 0j / j ^ j 1 ,m' z “c 1 j' ,m z 


(A.l) 
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with | ty \ rn ,, t^, | being the translation-addition coefficients, that can be found in [52, 53]. The magnetic 
field can be computed using the relation [54]: 


H = * 


V x E 


(A.2) 


Due to the symmetries of the problem, the interior and scattered field are given by [44]: 


V d vf (e) A (e) + r A (m) 

(A.3) 

V o-,f (e) A (e) +b-,t {m) A (m) 

j',rn' z 

(A.4) 


where {a^, bj>, Cj>, dj>} are functions of the multipolar order j; the ratio between the index of refraction of 
the resonator and the embedding medium n r ; and the size parameter of the problem, given by x = , 

with R the radius of the spherical resonator and A the wavelength of the emission of the radiating dipole. 
The coefficients a,, bj , Cj , dj are known as Mie coefficients when the dipole is outside of the resonator, but 
take a different form when the dipole is inside of it [23, 25]. The magnetic fields H lnt and H s are computed 
with eq.(A.2). 
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